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1.  INTRODUCTION 


For  an  airborne  antenna  mounted  on  the  fuselage  of  an  arbitrary' 
aircraft  near  the  extreme  top  or  bottom,  the  fuselage  shape  has  the 
dominant  effect  on  the  resulting  radiation  pattern.  The  wings  and  other 
flat  plate  structures  can  have  strong  effects  in  certain  sectors  of  the 
pattern  but  they  are  not  as  dominant  as  the  fuselage  especially  when 
the  complete  volumetric  pattern  is  concerned.  Previous  analyses  fl,2] 
for  the  airborne  antenna  performance  employed  a  profile  cylinder  model 
for  the  evaluation  of  the  radiation  pattern  in  and  near  the  elevation 
plane;  whereas,  a  cross-section  cylinder  model  was  used  for  the  pattern 
calculation  elsewhere.  Thus  for  pattern  calculations  off  the  principal 
planes,  a  two-cylinder  model  (21  was  employed  for  the  analysis.  Although 
this  two-cylinder  model  provides  a  useful  means  to  obtain  satisfactory 
results  for  the  volumetric  pattern  calculat ions,  there  are  some  problems 
encountered  in  the  nose  and  tail  sectors  of  the  pattern.  This  results 
from  the  fact  that  the  geodesics  (ray  paths)  used  in  the  twc-cyl  inc'er 
model  does  not  travel  on  the  true  three-dimensional  surface  of  the  fuse¬ 
lage.  As  a  consequence,  the  unit  normal  and  binormal  vectors  of  the 
geodesic  for  the  two  cylinders  do  not  vary  continuously  on  the  surface. 
This  leads  to  erroneous  results  and  discontinuities  in  the  calculated 
patterns'.  Due  to  the  above  mentioned  difficulties,  it  is  apparent  that 
an  improved  three-dimensional  model  is  needed  for  the  fuselage  structure. 
An  attempt  has  been  made  previously  to  model  the  fuselage  as  a  body  of 
revolution  [1],  and  good  results  were  obtained  for  an  antenna  radiating 
from  a  prolate  spheroid.  According  to  the  geometrical  theory  of  dif¬ 
fraction  (GTD)  analysis,  an  antenna  located  on  the  surface  launches  a 
ray  which  propagates  along  the  geodesic  path  on  the  surface  and  contin¬ 
uously  sheds  electromagnetic  energy  in  its  tangent  direction.  A  typical 
ray  picture  for  an  antenna  radiating  from  a  surface  of  revolution  is 
shown  in  Fig.  1,  in  which  a  ray  launches  from  S,  and  traverses  along 
the  geodesic  path  to  the  shedding  point  Q,  and  then  propagates  along 
the  tangent  direction  at  Q  in  the  observation  direction  r(e^,(|)  ).  Thus 
a  major  task  in  the  analysis  is  to  determine  the  unique^geooesic  path 
SQ  associated  with  a  given  desired  radiation  direction  rf©  ,4)^.).  The 
previous  numerical  technique  fl]  was  not  very  efficient  in'^that  a  tedious 
interpolation  procedure  was  employed  to  determine  the  unique  geodesic 
path  whose  tangent  points  in  the  desired  radiation  direction. 

In  this  report,  a  three-dimensional  model  for  the  fuselage  shape, 
namely,  a  body  of  revolution  model,  is  re-examined  and  an  efficient  method 
is  developed  to  find  the  unique  geodesic  corresponding  to  the  desired 
radiation  direction.  This  method  is  described  in  detail  in  the  following 
section. 

II.  GEODESICS  ON  A  BODY  OF  REVOLUTION 

Consider  a  body  of  revolution  with  the  parametric  representation 
of  its  surface  given  as  follows: 
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(1) 


X(",<:)  =  Rsin-cosifx  +  Rsin'-sin:y  +  Rcos'Z 


where  {R,e,<})  is  the  usual  spherical  coordinate  system,  and  R  is  a  func¬ 
tion  of  theta  i.e.,  R=R(o)-  The  surface  is  illustrated  in  Fig.  2. 

It  can  be  shown  that  tiie  arc  length  between  two  points  on  the  bodv 
of  revolution  can  be  expressed[33  as 

i  =  J  LR'^+  +  R^sin^fl((;‘ )•  ]^''^dr-  (2) 

iri  which  R'  and  are  the  derivatives  with  respect  to  the  variable  P. 
According  to  the  calculus  of  Variation[4] ,  ttie  arc  length  with  minir.uni 
distance  between  the  two  points  (the  geodesic)  should  satisfy  the  following 
constraint; 


R^sin^gi;' _ 

[R'^+R^+R^sin^e(4>‘ )‘  ]^ 


constant  =  A. 


(3) 


Equation  (3)  can  be  rearranged  and  one  obtains  the  following  expression 
for  the  geodesic  on  a  body  of  revolution: 


or 


4'  =  A 


(R' )^+R^ 


R  ^  i^o7^s7n^e-A^ 


1/2 


y  a 


,1/2 


(R')^+r^ 


?  ?  7  9  7” 

R  sin  e (R^sin^e-A^) 


dP 


(4a) 


(4b) 


In  the  above  equation  {q  ,4>  )  is  the  starting  location  of  the  geodesic, 
and  the  constant  A  is  reTatld  to  the  starting  direction  |  at  the  start¬ 
ing  location  through  Equation  (3). 

As  an  example,  let  us  consider  the  prolate  spheroid  v/ith  major  and 
minor  axes  "a"  and  "c"  as  shown  in  Fig.  3.  For  the  prolate  spheroid 
shown  in  Fig.  3,  it  is  easy  to  show  that 


R(e)  = 


_ ac _ 

2  2  2  2  ^ 

(a  cos  e+c  sin  e) 


(5) 
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and 


R' 


2  2 

-ac  sinecose(c  -a  ) 

,2  2  2.2 

(a  cos  6+c  sin  e) 


(6) 


By  specifying  e  ,<j!  and  4'' the  starting  position  on  the  prolate 
spheroid,  Equations  (4^,  fS)  and  (B)  can  be  employed  to  generate  a  family 
of  geodesic  curves  on  the  spheroid.  It  is  convenient  to  introduce  a  new 
variable  such  that 


4'  I  ^  =  -cottOg 


(7) 


Substituting  Equation  (7)  into  (3),  one  obtains 


A  = 


R^sin^p(-cotw^) 

[(R')^R^+R^sin^rCOt^o 


(8) 


Notice  that  for  w  equal  to  90°  or  270°,  the  geodesic  equation  reduces 
to  4=4  ,  i.e.,  thi  geodesic  is  the  generating  curve  of  the  body  of  revolu¬ 
tion,  the  meridians.  It  should  also  be  noted  that  Equation  (4)  does  not 
predict  the  geodesic  circle,  the  Equator,  e=90°  as  a  result  of  choosing 
Equation  (4)  as  the  representation  of  the  geodesic.  Now  let  us  examine 
some  numerically  computed  geodesic  paths.  In  Fig.  4  the  geodesic  curves 
for  a  starting  location  e  =90°,4  =0°  on  a  2x4  prolate  spheroid  are  plotted. 
For  a  given  staring  direction  w  ?  the  geodesic  path  is  specified  by  a  set 
of  points  on  the  spheroid  surface  with  e  and  4  being  the  position  parameters 
In  Figs.  5  and  6,  geodesic  curves  are  shown  for  two  different  starting 
locations,  (0  =60“, 4  =0°)  and  (e  =30°,4  =0°),  respectively.  Since  the  air¬ 
craft  fuselage  usually  resembles  a  long  cylindrical  structure,  Fig.  7  shows 
the  geodesic  curves  for  a  2xl00  prolate  spheroid.  It  is  interesting  to 
notice  that  the  geodesic  curves  for  the  2x4  and  2x100  spheroid  are  very 
similar  to  each  other.  Finally,  the  geodesic  for  a  sphere  is  shown  in  Fig.  8 

Now  that  the  geodesic  paths  are  known  for  the  surface  of  revolution, 
the  radiation  direction  along  the  geodesic  must  be  found  for  each  path. 
Recall  that  the  parametric  representation  of  the  geodesic  on  a  surface 
of  revolution  is 
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y  = 


Rsin6cos<ti  X 
Rsin6sin<p  ^ 
Rcose  z 


wi  th 


<^3  + 


R'^+R^ 


■lV2 


[V^Tn^e  { R^s  i  n^e -A^ ) 


de 


(8) 


(9) 


The  unit  tangent  vector  along  the  geodesic  is 


t  = 


dT  _  dy  di 
di  ”  de/de  ' 


It  can  be  shown  that  Equation  (10)  can  be  expressed  as  follows: 


(10) 


where 


^  _  cos(t>(R' sine+Rcose)-4>' siniiRsine 
■  0 


(11) 


(12) 


^z  = 


s in(ii(R'  si ne-t-Rcose )+())*  cos^Rsi ne 


R'  cos0-Rsine 


(13) 

(14) 


D  =  ^  =  sign(e-e^)j^^+R^+R^sin^e(^' )^  (15) 

Note  that  <f>'  is  given  by  Equation  (4a). 

It  is  convenient  to  express  the  tangent  direction  of  the  geodesic 
in  terms  of  the  spherical  coordinate  parameters  and  $  .  From  Equa¬ 
tions  (12),  (13)  and  (14),  the  tangent  direction  angles  6^  and  (fi  to 
the  geodesic  path  can  be  determined  as  follows:  ^ 


1? 


=  cos' 


-  tdn 


Numerical  results  for  ttie  tangent  direction  in  terms  of  e.  and 
4..  for  various  geodesic  curves  are  shown  in  Figs.  9,  10  and  11,  ^It  is 
interesting  to  observe  that  for  the  sphere  case  tlie  geodesic  tangent 
curves  converge  to  one  point,  for  the  circular  cylinder  case  tiiese  curves 
are  parallel  to  each  other,  and  finally  for  the  prelate  splieroid  case, 
the  tangent  curves  intersect  each  other  in  a  caustic  region.  Once  the 
unit  tangent  vector  is  deterniirK  <!,  the  other  relevant  parameter  in  the 
ray  analysis,  namely,  the  unit„ binprmal  vector  6  can  readily  be 
obtained  via  the  equation  b  =  t  x  n.  The  unit  vector  n  is  the  outward 
normal  vector  of  t!ie  surface. 

Determination  of  the  Geodesics  on  a  Prolate  Spheroid 

In  the  proceeding  section,  a  family  of  geodesic  curves  correspond¬ 
ing  to  a  particular  starting  location  on  the  surface  of  revolution  has 
been  determined.  The  associated  tangent  direction  curves  for  each  geo¬ 
desic  has  also  been  found.  According  to  the  geometrical  theory  of  dif¬ 
fraction,  energy  propagates  outward  from  the  source  along  the  surface 
geodesic  and  continuously  sheds  electromagnetic  energy  in  its  tangent 
direction.  Therefore,  the  next  step  is  to  determine  the  unique  shedding 
point  on  a  particular  geodesic  path  at  which  the  tangent  points  in  the 
desired  radiation  direction.  This  task  is  performed  in  this  section. 

Consider  now  an  antenna  is  located  on  a  surface  of  revolution  at 
S  as  shown  in  Fig.  1.  The  location  parameters  for  the  source  point  are 
The  desired  radiation  direction  (ej.,()>  )  is  denoted  with  r 
whfen  can  be  decomposed  into  its  rectangular^components  as  follows: 


TxX  +  r 


where 


Tx  =  sine|,cos<f>^ 


Ty  =  sine^sin4i^ 


=  cose^ 
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Recall  that  the  unit  tangent  vectors  for  the  geodesic  are  given  by 
Equations  (11)-(15).  Assuming  at  a  particular  point  on  the  sur¬ 

face  of  revolution,  the  tangent  direction  tU  matches  witH  the  radia¬ 
tion  direction  r,  then  from  Equations  (ll)-tl5)  and  Equations  (18)-(21) 
the  following  expressions  are  obtained: 


cositglR'  sineg+RcoseQ)-sin<QRsineg({.' 

“d7 


=  sine^cosit^ 
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<fy  -  COS 


-1 


cose„  R‘sinerv+Rcose„ 

_ _ _ 

sine^  R'coseg-Rsinep  ^  ^ 


(2B) 


In  Equations  (27)  and  (28),  there  are  three  unl.nowns  fp,  and  Aq 
and  one  needs  one  more  equation  to  solve  for  the  unknowns.  This  final  ^ 
equation  is  provided  by  the  geodesic  equation.  Equation  (4b),  derived  in 
Section  II,  namely. 


2  2 
R'^+R 


R  sin^e(Rsin^e-AQ^) 


(29) 


From  Equations  (28)  and  (29)  one  obtains  the  following  expression: 


^s  "  \ 


(30) 


Let  us  now  proceed  to  eliminate  Aq  from  Equation  (30).  From  Equation  (28) 
one  obtains  the  following  equation: 


cos(*  -<^n)  =  77 


cose. 


r  ^0'  sine. 


R'sineq+Rcoseq 
R ' coseq-RsinBq 


(31) 


Let  us  define 


cosa  =  COs(4)|,-(tiq), 


(32) 


then  from  Equation  (27),  Aq  can  be  written  as: 
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Aq  =  sign(eQ-e^)RsineQsine^/l-cos^ci 


(33) 


Note  that  cosa  is  related  to  Gq  via  Equation  (31)  and  Equation  (32), 

Finally,  Equations  (30)  through  (33)  can  be  employed  to  determine 
the  unique  geodesic  path  and  the  unique  diffraction  point  for  a  given 
source  location  on  the  surface  of  revolution  and  a  given  desired  radiation 
direction.  The  secant  method  is  suitable  for  solving  the  nonlinear 
equation.  Equation  (30).  The  initial  points  employed  for  the  secant 
method  are  and  c.  where  U  c'ven  by  Ecua'.icr 

A  simple  and  efficient  computer  program  has  been  developed  to  solve 
the  nonlinear  equation.  Equation  (30)  for  the  unknown  Gq.  The  secant 
method  with  the  initial  points  e  and  e^,  has  been  shown  to  be  an  extremely 
efficient  and  accurate  method  which  converges  in  a  few  iteration  steps. 

Once  the  is  determined,  Aq  and  fg  can  be  found  via  Equations  (33)  and 
the  'geodesic  equation.  Equation  (28;,  respectively.  Thus  for  a  given 
source  location  (65,41^),  the  unique  geodesic  path  with  a  starting  direc¬ 
tion  Aq,  and  the  unique  diffraction  point  (©qi^q)  on  the  surface  of  rev¬ 
olution  that  will  radiate  toward  the  desired^raaiation  direction  (e^,$|.) 
has  been  found  in  an  efficient  manner. 

Conclusion 

A  major  task  in  the  application  of  the  geometrical  theory  of 
diffraction  (GTO)  to  the  problem  of  electromagnetic  radiation  from  a 
general  convex  surface  is  to  determine  the  unique  geodesic  path  that 
starting  from  the  source  location  traverses  the  surface  and  has  a  tangent 
at  some  point  along  the  path  which  points  in  the  desired  radiation  direc¬ 
tion.  In  this  report  a  numerically  efficient  and  accurate  scheme  has 
been  developed  for  solving  this  problem  in  the  case  of  a  general  convex 
surface  of  revolution.  The  surface  of  revolution  is  of  interest  in  that 
it  provides  an  analytic  model  for  the  aircraft  fuselage  structure.  A 
computer  program  is  also  developed  to  solve  the  governing  nonlinear 
equation  using  the  secant  method  for  iteration.  This  method  is  proved 
to  be  simple  and  converges  to  the  accurate  results  in  just  a  few  steps. 
Numerical  results  for  the  family  of  geodesic  curves  are  also  presented 
for  the  cases  of  sphere,  prolate  spheroid  and  cylinder. 
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